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ABSTRACT 

We present the first multi-wavelength follow-up observations of two candidate gravitational-wave 
(GW) transient events recorded by LIGO and Virgo in their 2009-2010 science run. The events were 
selected with low latency by the network of GW detectors (within less than 10 minutes) and their 
candidate sky locations were observed by the Swift observatory (within 12 hours). Image transient 
detection was used to analyze the collected electromagnetic data, which were found to be consistent 
with background. Off-line analysis of the GW data alone has also established that the selected GW 
events show no evidence of an astrophysical origin; one of them is consistent with background and the 
other one was a test, part of a "blind injection challenge". With this work we demonstrate the feasibility 
of rapid follow-ups of GW transients and establish the sensitivity improvement joint electromagnetic 
and GW observations could bring. This is a first step toward an electromagnetic follow-up program in 
the regime of routine detections with the advanced GW instruments expected within this decade. In 
that regime multi-wavelength observations will play a significant role in completing the astrophysical 
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identification of GW sources. We present the methods and results from this first combined analysis 
and discuss its implications in terms of sensitivity for the present and future instruments. 
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1. INTRODUCTION 

Some of the key questions in the pursuit of sources 
of transient gravita tional radiat ion detectable by LI GO 
(| Abbott et al.ll2009D and Virgo (jAccadia et al.ll2012l ) re- 
late to their electromag netic (EM) signature s and our 
ability to observe them (lAbadie et al.l 12012 al) . In sev- 
eral of these sources, like core collapse supernovae and 
neutron star-neutron star or neutron star-black hole 
mergers, energetics suggest gravitational waves (GWs) 
are likely to be accompanied by EM emission across 
the spectrum and over time scales ranging from sec- 
onds t o days (iFrver et al l 120021: iPiran Il2004t iMeszaros I 
[200l INakar I 120071 : ICorsi fc Meszarosl \mm . Multi- 
wavelength EM observations of such events have al- 
ready set the paradigm for improved constraints on 
source astrophysics set jointly rather than separately. 
Prompt outbursts, as well as afterglows in X-ray, optical 
and radio associated with Gamma- Ray Bursts (GRBs) 
and supernovae have shed light o n the progenitors and 
the astrophysics o f these systems ([Kulkarni et a.1.1 Il998t 
Bloom et al.l 11999b [Matheson et al.l 120031: [Berger et al. I 
20051: [Gehrels et al.ll2005HSoderberg et al.ll2005l iBerger l 
200a IZhang fc MacFadvenl 120091: iSoderberg et aJJl201Ct 
Berger 1 1201 W ~ The benefit of coupling them to GW 
observations will be tremendous as this will bring 
the GW observations into astrophysical and cosmo- 
logical context. Besides incr easing detection confi- 
dence (jKochanek fc Piranl fl993). multi- wavelength ob- 
servations may improve source localization down to 
the arcsecond level, leading to identification of the 
host galaxy and meas u rement of the r edshift dSchutz I 
llMlSyT^estrel i| l200l iStubbs I 120081 iPhinnevI 120091: 
Stamatikos et al.ll2009t iBloom et al.ll2009l : iMetzger et all 
2010HMetzger fc Berger II2012I ). Conversely, the absence 
of any EM signature for an otherwise confident transient 
detection with the GW detectors alone will provide con- 
straints on emission mechanisms, progenitors and ener- 
getics. 

LIG O (| Abbott et al.ll2009f ). and Virgo (jAccadia et al.l 
[2012J) form a network of interferometric detectors aim- 
ing to make the first direct observations of GWs. In 
their 2009-2010 data-taking period this network con- 
sisted of three interferometers: LIGO-Hanford in Wash- 
ington State in the USA, LIGO-Livingston in Louisiana 
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in the USA, and Virgo in Italy. During these science runs 
of the instruments we implemented a first program that 
could allow prompt EM follow-up of candidate GW tran- 
sients. Starting with GW data, lo w-latency searches for 
compact binary star coalescences (jAbadie et al.ll2011a| ) 
and un-modeled GW transients were performed. This 
allowed the prompt identification and sky localization of 
GW candidates, which for the first time were passed on 
to ground-based telescopes and Sw ift in order to be fol- 
lowed up. In an earlier publication {Abadic ct al.ll2012aT) 
we presented the details of the implementation and test- 
ing of this low latency search and its coupling to EM 
astronomy. In this paper we report on the analysis of 
the Swift data we collected as part of this program. 

With the start of the 2009-2010 LIGO- Virgo runs we 
established a Target-of-Opportunity (ToO) program with 
Swift in order to search for possible afterglow in X-ray, 
ultraviolet and optical wavelengths of a small number 
of GW transient candidates. The sensitivity of the GW 
detectors at the time made the chance of a detection 
small, but non- negligible. For the case of binary sources 
with at least one neutron star in the system, the rate 
of detectable merger events was predicted to be in the 
range of 2.7xl0~ 4 to 0.3 events per year (jAbadie et al.l 
l2010a|) . This ToO program was exercised twice; neither 
time led to detection of an EM counterpart to a GW 
transient. Nonetheless, the program addressed imple- 
mentation questions and established the first joint obser- 
vation and coordinated data analysis by the LIGO- Virgo 
network and Swift satellite. 

In this paper we present the results from the EM 
follow-up program involving only Swift — results from 
the sister program involving the follow-up of GW candi- 
dates by ground-based optical and radio telescopes will 
be the subject of a forthcoming publication. The paper 
is organized as follows. In Section 2 we describe the Swift 
observatory. We then review in Section 3 the procedure 
for targetin g EM follow-up of GW events, as described 
in detail in lAbadie et al.l (|2012al) . The Swift observa- 
tions and analysis of data are described in Section 4. In 
Section 5 we present our formalism for combining results 
from the joint GW-EM search, including the simulation 
work we performed. We conclude in Section 6 with a 
discussion and outlook for this kind of joint search. 

2. THE SWIFT OBSERVATORY 

Th e Swift Gamma-Ray Burst Mission () Gehrels et al.l 
I2004D . developed and launched under NASA's Medium 
Explorer Program, is unique for its broad wave- 
length sensitivity and rapid response. Three tele- 
scopes are co-aligned. The Burst Alert Telescope 
(BAT) (|Barthelmv et al.1 12005) is a broad field of view 
(FOV) coded-aperture instrument with a CdZnTe detec- 
tor plane, designed to search for transient events such as 
GRBs. Approximately 100 GRBs are discovered by BAT 
each year. In response to a BAT trigger, the spacecraft 
performs an autonomous slew to point the two narrow 
field of view instru ments. The Ultra- V iolet and Optical 
Telescope (UVOT) (|Roming et al.ll2005l ) performs follow- 
up observations of GRBs in the 170 to 600 nm band, 
with a 0.28° x 0.28° field of view. The X-Ray Telescope 
(XRT) (jBurrows et al.lf2005h has an effective area peak- 
ing at 110 cm 2 (at 1.5 keV), a 0.4°x0.4° FOV, and an 
energy bandpass of 0.3-10 keV. While spacecraft slew- 
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ing can be initiated autonomously, achieving repointing 
within approximately 2 minutes, the process can also be 
initiated from the ground. The Swift Mission Operations 
Center (MOC), located near the University Park cam- 
pus of the Pennsylvania State University, receives over 
1000 requests for ToOs from the astronomical commu- 
nity each year. Response time to a ToO depends on 
scientific priority and urgency. Due to communication 
limitations and human-in-the-loop commanding, obser- 
vation of the highest priority ToOs is typically achieved 
within 4 hours, although frequently the response time is 
less than 1 hour. In this way, a ground-based observa- 
tory such as LIGO- Virgo can provide the trigger for a 
transient event, with Swift providing prompt follow-up 
observations in the X-ray, UV, and optical bands. 

3. SELECTION OF GRAVITATIONAL- WAVE TRANSIENTS 

The LIGO- Virgo EM follow-up program took place 
during two observing periods spanning a bit over two 
months from 2009 December 17 to 2010 January 8 and 
from 2010 September 2 to 2010 October 20. The de- 
tails of the implementation and testing of the end-to-end 
search on the GW e nd have been described elsewhere 
(jAbadie et alJl2012af ). Here we will summarize the key 
features and expand on aspects that were specific to the 
joint program with Swift. 

3.1. Event selection 

GW candidate events were selected for follow-up based 
primarily on their False Alarm Rate (FAR), the rate at 
which an event of equal or greater significance is expected 
to occur in the absence of a true signal. For this run, 
in order to trigger a Swift follow-up, we set a nominal 
threshold on the candidate events' FAR at no more than 
one event per 35 days of triple coincident running (i.e., 
the GW instrument configuration when all three detec- 
tors were acquiring sensitive data). 

To avoid trying to image unusually poorly localized 
events, an additional constraint was placed on candidates 
that 20% or more of the a posteriori (see Section 13.21) 
weighted probability in the event skymap must be cov- 
ered by up to five 0.4° x 0.4° tiles selected for follow-up 
as Swift fields. The number of fields picked for the Swift 
follow-up program reflected a compromise between the 
need to cover as large an area in the sky as possible 
and the requirement to be minimally disruptive of Swift's 
other science targets. If five fields were targeted by Swift 
for each LIGO- Virgo candidate, the correct location of 
GW events of high significan ce could be imaged w ith a 
probability greater than 50% (jAbadie et al.ll2012at ). Ad- 
ditionally, manual and automated checks were performed 
on data quality in each interferometer before sending any 
alerts, eliminating events which had an obvious problem 
associated with them. For events which passed all data 
selection criteria, observation requests were sent to Swift 
through web-based ToO submissiorQ. 

During the course of the EM follow-up program, GW 
candidate transient events were selected by transient- 
finding algorithms for un-mo deled bursts as well a s com- 
pact binary star coalescences (jAbadie et a l. 2012a). Two 
such events met all criteria and were submitted for follow- 
up with Swift. These were identified by one of the generic 

1 https:/ /www. swift. psu.edu/secure/toop/too_request. htm 



transient-finding alg orithms, called coherent WaveBurst 
(jKlimenko et al.ll2011l) . This algorithm uses wavelet de- 
composition to search for GWs without relying on spe- 
cific waveform models. The "January" event occurred at 
8:46 UTC on 2010 January 7. This was an event close 
to the end of the first LIGO- Virgo observation period, 
during which we lowered the nominal thresholds to 1 
event per day for the FAR and 10% for the weighted 
probability of the event skymap. The thresholds were 
adjusted as the first LIGO- Virgo observation period was 
approaching an end in order to exercise the follow-up 
process at least once. The "September" event, which oc- 
curred at 6:42 UTC on 2010 September 16, passed all 
nominal criteria for follow-ups. This event was later (in 
March 2011) revealed to be a "blind injection" artificially 
inserted into the interferometers as a te st of o ur detec- 
tion procedures (jAbadie et all I2011bl lc_ I2012H ). While 
both of these events were ultimately tests of the system 
rather than plausible candidates, they demonstrate the 
viability of performing rapid follow-ups of potential GW 
signals using Swift. 

3.2. Position errors and tiling 

The typical uncertainty in sky location of a GW signal 
is large (typically tens of square degrees) relative to the 
FOV of Swift's XRT and UVOT instruments. This may 
be in part addressed by imposing the requirement for the 
reconstructed sky location to overlap with nearby galax- 
ies. We have used information f rom the Gravitatio nal 
Wave Galaxy Catalog (GWGC) (IWhite et al.ll20ll in 
order to fold into the tiling algorithm location, extent 
and blue luminosity of known galaxies within a distance 
less than 50 Mpc, since the GW interferometer network 
would not be likely to detect neutron star binary coales- 
cences beyond this distance. In this way, locations with a 
known galaxy from this catalog were given greater weight 
than those without galaxies present. 

For the earlier "winter" run in which the January event 
occurred, skymap tiles with galaxies or globular clusters 
at a distance less than 50 Mpc had their estimated prob- 
ability increased by a factor of 3, whereas for the later 
"autumn" run the estimated relative probability of each 
tile in the skymap was calculated according to: 



where Mi is the blue lig ht luminosity of a galaxy (a proxy 
for star formation rate dPhinnev lll99lHKopparapu et al.l 
l2008HAbadie et alll2010al) L L is the likelihood from only 
the GW sk ymap and Dj is the dist ance of the galaxy 
from Earth (jNuttall fc Sutton |[2010h . The index i sums 
over each galaxy associated with the skymap tile. Only 
~8% of tiles in a typical skymap were associated with 
one or more galaxies, whil e P was set to zero for tiles 
with no associated galaxy (jAbadie et al.ll2012af ). After 
this summation was performed, the resulting likelihoods 
were renormalized to a probability of unity over the en- 
tire skymap. In this way we constructed the a posteri- 
ori probability skymaps that we used for prioritizing the 
Swift observations. 

For both events selected for follow-up, several tiles 
were chosen and passed on for imaging as Swift fields 
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Fig. 1. — Probability skymaps of the January event (top) and 
September event (bottom) as determined by the coherent Wave- 
Burst algorithm before galaxy weighting. The asterisks mark the 
approximate locations of the locations selected for follow-up with 
Swift. Color scale shows fraction of total (raw) probability con- 
tained within that particular skymap tile. 

according to the procedures outlined above. The prob- 
ability skymaps produced by coherent WaveBurst with- 
out galaxy weighting and containing the top 1000 tiles 
are shown in Figure [1] The maps show large spreads in 
probability on the sky, the first due to their relatively low 
significance, the second due to worse sensitivity in Virgo 
relative to the LIGO interferometers. Simulations show 
that, despite the relatively small size of Swift fields with 
respect to the extent of our probability skymaps, sources 
originating in or close to nearby galaxies ca n be correctly 
localize d with reasonable success rates (see lAbadie et al.l 
(|2012aD and discussion in Section 5). The regions con- 
taining the fields imaged by Swift are shown with an as- 
terisk on each skymap, with precise coordinates of each 
field given in Tables [THU 

4. OBSERVATIONS WITH SWIFT 

The ToO requests were submitted manually via the 
web within 2 hours from the collection of GW data. The 
initial two requested pointings for each event were sub- 
mitted as "high priority" (but not the highest); this re- 
flected a compromise between capturing the early light 
and being minimally disruptive to Swift operations. 
Swift observations with the XRT and UVOT instruments 
of the specified targets were performed in the following 



~ 12 hours. Rapid follow-up could significantly increase 
the likelihood of a detection, as several types of potential 
GW/EM sources (e.g. short GRBs) have afterglows that 
fade below detectability on a timescale of ~ hours. The 
highest probability region of the January event skymap 
was observed by Swift in 5 overlapping fields (Figure [5]) . 
Each of these were observed twice that month. For the 
September event two disjoint Swift fields were observed 
(Figure [2]). Note that the fields shown are explicitly 
for the XRT analysis, whereas the UVOT has a slightly 
smaller FOV. Each September event field was observed 
on two days, one in September 2010 and the other in 
December 2010. Details of the observations are given in 
Tables Q3- II 

4.1. X-ray results 

We analyzed the Swift-XRT data with custom scripts 
whic h use the so ft ware described in Evans et al. (2007) 
and lEvans et all (|2009f ). For each field we combined 
the observations and produced a single image and ex- 
posure map. We then used the source detection and 
point sprea d function ( PSF) fitting code described in 
iGoad et al.1 (|2007t ) and lEvans et al.l (|2009D to identify 
and localize sources in the field. This method uses a 
sliding-cell detection algorithm, with a fixed cell size of 
21x21 pixels (49.6"x49.6", which encloses 93% of the 
PSF). An initial run with a detection significance thresh- 
old of 3 (i.e. S/as is at least 3, where S is the estimated 
number of net source counts in the cell, and as is the un- 
certainty in this value, determined using the background 
estimated from an annular box of size 51 pixels) revealed 
no sources in any of the observations, however for faint 
sources this box size is sub-optimal. We performed a 
second source search with a reduced detection thresh- 
old of 1.5cr. We measured the mean background level 
in each field and, using a circular source region of ra- 
dius 10 pixels centered on the position of each of these 
"reduc ed-threshold " detec tions, we applied the Bayesian 
test of iKraft et al.l (|1991l) . only accepting sources which 
this test determined to be detected with at least 99.7% 
confidence (i.e. 3cr detections). While this method makes 
us more sensitive to faint sources, it also increases like- 
lihood of getting false positives due to background inho- 
mogeneities, thus these detections should be treated with 
caution. A detection system optimized for faint sources 
is under development. The positions of the reduced- 
threshold detections are given in Tables [3] |U 

Each of these reduced-threshold detections had very 
few photons, therefore it was not useful to perform a 
detailed spectral analysis. We determined the mean 
co unt-rate of e a ch de tection using the Bayesian method 
of IKraft et al.l ()199lD and then applied corrections for 
PSF losses and ins trumental effects, f ollow ing the pro- 
cesses described by lEvans et al.1 (|2007r i and lEvans et al.1 
( 2009]). We used PIMMS0 to determine a count rate-to- 
flux conversion factor, assuming an absorbed power-law 
spectrum with a photon index of 1.7. For the January 
event the assumed absorbing column was 4 x 10 20 cm -2 , 
giving a 0.3-10 keV conversion factor of 4.2 xl0~ n erg 
cm~ 2 ct _1 ; for the September event the absorbing col- 
umn was 3 x 10 21 cm -2 , giving a 0.3-10 keV conversion 
factor of 5.1 x 10~ u erg cm~ 2 ct _1 . A change of 0.2 in the 

2 http:/ /heasarc. nasa.gov/Tools/w3pimms. html 
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Fig. 2. — Swift-XRT image of the vicinity of the January event 
(top) and September event (bottom). The large circles show the 
fields of view of the various pointings. The grayscale denotes rel- 
ative length of exposure at each location with darker colors indi- 
cating longer ones. The horizontal axis denotes right ascension 
(RA) in hours, minutes and seconds and the vertical axis denotes 
declination (dec) in degrees. 

assumed index of the power-law spectrum would result 
in a ~10% change in these conversion factors. The ab- 
sorbing columns were taken as the Galact ic values in the 
direct ion of the event, determined from IKalberla et al. 
(2005) assuming the abundances of lAnders fc Grevesse 
(1989|). 

Given the sky area covered by the Swift observations 
(0.126 square degrees per XRT field after instrumental 
corrections are applied) we expect to find a number of 
serendipitous sources in the XRT data. To quantify the 
likelihood of serendipitou s source detection, w e used the 
2XMMi-DR3 catalogue (| Watson et all I2009D . which is 
substantially dominated by serendipitous X-ray sources 
found by the XMM-Newton observatory (|Jansen et al.l 



120011) . 

We selected from this catalogue all unique good sources 
(i.e. with a q u ality flag of 0; for the flag definition see 
IWatson et all (|2009D ). For consistency with 2XMMi- 
DR3 we convert the 0.3-10 keV XRT count-rates into 
0.2-12 keV fluxes, using an absorbed power-law spec- 
trum with a column of TVjj = 3 x 10 20 cm~ 2 and a power- 
law photon index V — 1.7, as used for 2XMIvfl. For each 
XRT source we then counted the number of 2XMMi-DR3 
sources with values of the 0.2-12 keV flux — as meas ured 
with the EPIC MOS-1 camera ([Turner et al.ll200l on- 
board XMM-Newton — at least as bright as the source 
in question. This was then scaled by the ratio of the 
Swift-XRT instantaneous field of view to the 2XMMi- 
DR3 unique sky coverage area (504 square degrees). This 
yielded, for each Swift source, the number of serendipi- 
tous sources (N s ) of at least that brightness expected in 
a single Swift field of view. These values, along with the 
count rates, are given in Tables [SHU For the January 
event we expected a total of 7.5 serendipitous sources, 
compared to our 8 detections; for the September event 
we expected 17 serendipitous sources, and had 12 de- 
tections. The number of serendipitous sources expected 
will not be a strong function of Galactic latitude at these 
latitudes (—27.6° for the January event and —6.1° for 
the September event) and flux levels; Galactic source sky 
densities are comparable to extragalactic values at a flux 
of around 5 x 10 ~ 14 erg cm" 2 s -1 o nly at much lower 
latitudes, see e.g. IMotch et ail (|2010D . 

In preparing for the event collection and Swift follow- 
up, we constructed a utility catalog that would allow 
us to look up the X-ray history of each field selected 
for observation in order to help determine if any of the 
observed sources coincide with known steady-state or 
variable sources. This X-ray source catalog, which we 
call XGWC0, combines public data from the HEASARC 
Master X-ray Catalog and the GWGC. Upon exami- 
nation of XGWC/HEASARC for known X-ray sources 
within the observed fields we found two X-ray sources 
(one of which was observed twice by Swift) in the cata- 
log which were within the fields for the January event. 
Both were of unknown type. These catalog sources are 
likely to be associated with the reduced-threshold sources 
#3 and #4 reported in Table [31 The fields that were ob- 
served for the September event contained only one XG- 
WC/HEASARC source of unknown type. This is likely 
to be associated with the reduced-threshold source #8 
reported in Table Q] 

We also performed a variability test on the observed 
sources, since an EM counterpart to the LIGO- Virgo 
event may be expected to be fading. To do this we pro- 
duced a light curve of each source, creating one bin per 
Swift observation IE0, and determining the count-rate 
using the Bayesian method (Kraft et al. 1991). The sig- 
nificance of any variability was then found simply by di- 
viding the difference in the two bins by the errors in the 
two bins, added in quadrature. This is given in Tables[3]- 

3 http:/ /xmmssc-www. star. le.ac.uk/Catalogue/2XMM/UscrGuide 
_xmmcat.html#TabECFs 

4 http:/ /aquarius. cltc.hu/XGWC/indcx. html 

5 http:/ /heasarc. gsfc.nasa.gov/W3Browse/all/xray.html 

6 For the September event we excluded the second observation 
of each source, since this was very short. 
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31 two sources show moderate (between 2 and 2.5a) evi- 
dence for variability (January Source #1 and September 
Source #6). 

In addition to the observations reported above, a fur- 
ther follow-up observation was performed for the Jan- 
uary trigger on 2010 Jan 30. This observation, obs 
ID 00031592001, contained 5.4 ks of data, centered at 
RA=88.73°, dec=-40.96°, i.e. overlapping heavily with 
field 00031577. Including the data from this field causes 
Sources #1 and #4 to be registered as above-threshold 
detections. For Source #4, the count-rate of the final 
observation was almost identical to the first; inclusion of 
this result changed the variability significance from 0.75 
to 0.04cr. For Source #1, the addition of these data alters 
the variability only slightly, from 2.05 to 1.98cr. 

4.2. Optical and UV results 

The LIGO- Virgo target fields were observed by the 
Swift UVOT with the clear white filter (160 - 800 nm) 
or, if ground analysis indicated that the field contained a 
star too bright for observation in white, the broadband u 
filter (centered at 346.5 nm with full width at half max- 
imum of 78.5 nm) or narrow NUV uvm2 filter (centered 
at 224.6 nm with full width at half maximum of 49.8 nm). 

We first attempted to identify any UVOT counter- 
parts to the X-ray reduced-threshold sources identified 
in the XRT images, photometering any source within 
5" of the XRT position using the standard UVOT- 
SOU RCE code and calibr ations from lPoole et all (|2008l) 
and IBreeveld et all (|2011l ). The results are listed in Ta- 
bles [SHSl For the January event, XRT reduced-threshold 
detection 7 fell off the smaller UVOT field and was not 
measured. UVOT counterpart sources to XRT reduced- 
threshold detections 1, 2, 3 and 4 were found in at least 
one epoch. However, they had corresponding sources in 
the Digital Sky Survey (DSS) and showed no photomet- 
ric variation beyond the l-2cr level. No counterparts to 
XRT reduced-threshold detections 5, 6 and 8 were found 
in UVOT. 

For the September event, XRT reduced-threshold de- 
tections 1, 7 and 8 fell outside of the UVOT field. No 
counterparts to XRT reduced-threshold detections 2, 4, 
11 and 12 were found by UVOT. Reduced-threshold de- 
tections 3 and 10 were also found by UVOT and cor- 
respond to DSS sources. XRT reduced-threshold detec- 
tions 5, 6 and 9 correspond to very marginal detections in 
a crowded field and are either very faint objects or spu- 
rious detections. The apparent dimming of the UVOT 
counterpart to source 5 is spurious as this source was on 
the edge of one UVOT image and only a partial flux could 
be measured. In short, no XRT reduced-threshold detec- 
tion corresponds to either an optical transient or a vari- 
able source for either the January or September event. 

We next performed a blind search for variable targets 
within the field. In this instance, we photometere d the 
entire field using the DAOPHOT (IStetson I U9871 PSF 
photometry program. We photometered all data us- 
ing a quadratically variable PSF and found photometry 
and image subtraction to be excellent, provided the PSF 
was based on stars not near the coincidence-loss limit 
of the data. Raw photometry was corrected for coin- 
cidence loss, exposure time and large scale sensitivity. 
The number of detections in each field ranged from 250 
to 2800 depending on the Galactic latitude of the field 



and the exposure time. Of the nearly 6800 sources that 
were detected in the UVOT data, approximately 5200 
are well-measured point-like sources that are not near the 
coincidence- loss limit off the data. Of these, 11 are not 
near chip edges or bright stars but show very significant 
variability (> 7 a) between the two epochs. Their coordi- 
nates are shown in Table [7] Such significant variability is 
likely to reflect real phenomena, such as active galaxies or 
variable stars, which are expected to be detected in any 
deep field. However, none of these variable sources cor- 
responded to any X-ray detections and none show vari- 
ability beyond a few 0.1 magnitudes, which would be 
consistent with a tentative classification as normal vari- 
able stars or active galaxies. Further monitoring would 
be needed to determine their nature. 

4.3. Summary of EM findings 

The XRT analysis produced 20 reduced-threshold de- 
tections and the UVOT analysis identified nearly 6800 
sources in the follow-ups of the two events. However, all 
observations in X-ray, UV and optical bands were consis- 
tent with expectations for serendipitous sources and no 
single source displayed significant variability in the XRT 
or UVOT analyses. 

5. COMBINED GW-EM RESULTS 

Information from the EM observations associated with 
a GW candidate event will have to be combined with GW 
data in order to establish key quantities associated with a 
combined GW-EM search for transients, like event signif- 
icance or astrophysical reach. In this section we present 
a formalism for combining results from such joint obser- 
vations. For the purpose of its validation we performed 
simulations using possible models of GW and EM sig- 
nals. In the process we verified the overall search pro- 
cedure and estimated the increase in sensitivity of the 
search which resulted from the Swift follow-up observa- 
tions. Although this is presented within the context of 
the LIGO- Virgo-Swift search, it can be straightforwardly 
extended to other joint searches as well. 

In the way the search was conducted, a typical joint 
candidate event is characterized by the strength of the 
GW signal expressed in terms of the coherent network 
amplitude 77 which is the detection statistic itself for 
coherent WaveBurst. This quantity is generally pro- 
portional t o the signal-to-noise r atio and is described 
in deta il in iKlimenko et al.l ([201 1[ ) and in lAbadie et al.l 
(|2012bf) . Joint candidate events are also characterized 
by the measured X-ray flux, S, and the sky location of 
the X-ray source, SI = [RA,dec]. In addition to these, a 
skymap, i.e. the probability distribution for the location 
of a potential GW source, p m (0), is produced and used 
to select fields for imaging. We define the joint detec- 
tion statistic for an event as the logarithm of the joint 
likelihood ratio, sometimes also referred to as the Bayes 
factor, given by 

Ajoint 

(77, S, SI) = A GW ( V )A EM (S)A C or(n). (2) 

In the above equation Aqw — p(v I signal) /p(rj | noise) 
is the likelihood ratio for a GW candidate, measuring 
its significance. Aem(<5) = p$ 1 (S) is the inverse of the 
probability density of observing an accidental, serendip- 
itous X-ray source which is not correlated with the GW 
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TABLE 1 

Details of the Swift- XRT follow-up observations of the January event. 



Swift Observation ID 




Date Start 
CUT) 






Date End 
(UT) 






Exposure 
(ks) 


.Pointing Direction 
J2000 


00031575001 


2010 


Jan 07 at 


13 


;04 


2011 


Jan 07 at 


11 


29 


1.5 


89.07", 


-40.96° 


00031575002 


2010 


Jan 11 at 


19 


11 


2010 


Jan 11 at 


19: 


38 


1.6 


89.00°, 


-40.92° 


00031576001 


2010 


Jan 07 at 


14 


:33 


2010 


Jan 07 at 


15: 


55 


1.7 


89.14°, 


-40.80° 


00031576002 


2010 


Jan 11 at 


15 


:59 


2010 


Jan 11 at 


16: 


18 


1.1 


89.13°, 


-40.77° 


00031577001 


2010 


Jan 07 at 


15 


:57 


2010 


Jan 07 at 


17 


20 


1.8 


88.60°, 


-40.81° 


00031577002 


2010 


Jan 11 at 


17 


:35 


2010 


Jan 11 at 


18: 


02 


1.6 


88.60°, 


-40.76° 


00031578001 


2010 


Jan 07 at 


17 


23 


2010 


Jan 07 at 


18 


49 


2.0 


88.07°, 


-40.79° 


00031578002 


2010 


Jan 11 at 


13 


:18 


2010 


Jan 11 at 


15: 


01 


1.4 


88.06°, 


-40.78° 


00031579001 


2010 


Jan 07 at 


20 


:26 


2010 


Jan 07 at 


20 


59 


2.0 


88.71°. 


-41.16° 


00031579002 


2010 


Jan 11 at 


18 


:03 


2010 


Jan 11 at 


19: 


51 


l.i 


88.77°, 


-41.17° 



TABLE 2 

Details of the Swift-XRT follow-up observations of the September event. 
Swift Observation Date Start Date End Exposure Pointing Direction 



(UT) (UT) (ks) J2000 



00031825001 


2010 


Sep 


16 


at 


18:10 


2010 


Sep 


16 


at 


18:35 


1.5 


110.98", 


-27.53" 


00031825002 


2010 


Dec 


30 


at 


03:17 


2010 


Dec 


30 


at 


03:19 


0.09 


110.98°, 


-27.54° 


00031825003 


2010 


Dec 


30 


at 


00:06 


2010 


Dec 


30 


at 


03:34 


2.3 


110.98°, 


-27.53° 


00031826001 


2010 


Sep 


16 


at 


19 36 


2010 


Sep 


16 


at 


20:09 


2.0 


109.86°, 


-27.57° 


00031826002 


2010 


Dec 


29 


at 


03:37 


2010 


Dec 


29 


at 


04:50 


0.14 


109.85°, 


-27.54° 


00031826003 


2010 


Dec 


29 


at 


03:38 


2010 


Dec 


29 


at 


05:09 


1.9 


109.86°, 


-27.58° 



TABLE 3 

The reduced-threshold detections in the X-ray data for the January event. 



Source # 


Right ascension (RA) 


Declination (dec) 


Position Error 


Count Rate 


N s 1 


Variability a 




(J2000) 


(J2000) 


(" 90% conf.) 


(0.3-10 kcV, ks" 1 ) 




significance (<r) 


1 


05h 55m 1.00s 


-40°58' 00.8" 


4.5 




0.9 


2.05 


2 


05h 57m 4.80s 


-40° 54' 45.4" 


4.3 




0.9 


0.26 


3 


05h 54m 12.72s 


_40°44' 05.8" 


1.3 


4.611:1 


1.3 


0.45 


4 


05h 54m 59.29s 


-40°54' 19.6" 


4.5 


q 9+1.3 
•J' z -1.0 


2. 1 


0.75 


5 


05h 51m 57.66s 


-40°46' 10.9" 


5.6 


2.8+lf 


2.9 


1.10 


6 


05h 51m 41.12s 


_40°44' 46.4" 


5.5 


1 4+ 1 - 1 


7.5 


0.74 


7 


05h 52m 6.29s 


-40°59' 14.3" 


6.5 


o q+1-2 
z '°-0.8 


3.9 


0.91 


8 


05h 52m 55.88s 


_40°46' 14.9" 


5.2 


2 9+ 1 - 7 
z.y_ 1 2 


2.8 


2.00 



a N s is the number of 2XMMi-DR3 sources which are at least as bright as the Swift source, which are expected in a single Swift field. Sec 
text for details. 

b The significance of the difference in count-rate between the difference epoch observations. 



signal. The remaining term A cor (f2) = p m (f2) is the 
probability for a GW source to be in the location of the 
X-ray source, which measures positional correlation be- 
tween GW and EM signals. Detection statistics based 
on the likelihood ratio or the Bayes factor construction 
have been previously sugge sted and used in the context of 
searches for GW bur sts (jClark et al]|2007l : iCannon ||200H; 
lAbadie et al.ll2011df) and GWs from compact binar y co- 
alescence (jAbadie et afll2010bl : [Biswas et al.l I2012D . In 
deriving Equation ((2|) (see Appendix [X] for details) we 
assumed that the dominant background in the EM sector 
is serendipitous X-ray sources that happen to be within 
the observed fields by chance. We neglected contribu- 



tions to this from possible spurious sources - due to, 
e.g., instrumental artifacts - in the XRT analysis. In 
order to check if this assumption is .jus tified we esti- 
mate d from the 2XMMi-DR3 catalogue (jWatson et al.l 
I2009D that one expects to find at least one serendipitous 
source with flux equal or greater than 5.4 xl0~ 13 erg 
cm~ 2 s _1 within five Swift fields. Visual inspection of 
the XRT data would identify artifacts at or above this 
flux level, and they would be excluded from the analysis. 
No such artifacts were found in the fields analyzed in this 
paper. We also assumed that inhomogeneities in the dis- 
tribution of serendipitous sources over the sky are small. 
As already mentioned earlier, for every LIGO- Virgo can- 
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TABLE 4 

The reduced-threshold detections in the X-ray data for the September event. 



bource # 


Right ascension (RAJ 


Declination (dec) 


Position Error 


n i — n — i 

Count Hate 


— KT 


t7 — 1 — rn 

Variability 




( 19000*1 


(J2000) 


^ yu/o coin. j 


1 \J . O J- \J ivc V • ivS I 






1 


07h 23m 22.99s 


-27° 26' 


10.1" 


4.4 


2 8 +u ^ 


2.9 


1.47 


2 


07h 23m 22.34s 


-27° 33' 


09.5" 


4.4 


2 3+ii 

z -°-0.7 


3.9 


1.09 


3 


07h 23m 34.43s 


-27° 23' 


32.4" 


5.4 


— U.o 


3.7 


1.47 


4 


07h 24m 34.95s 


-27°31' 


31.1" 


6.1 


1 8+ 1 ' 2 
1 - 8 -0.7 


5.5 


1.01 


5 


07h 23m 53.50s 


-27° 23' 


06.5" 


4.4 


6+ ' 3 
U.D_ 2 


17 


1.30 


6 


07h 24m 27.89s 


-27° 35' 


40.8" 


6.5 


2 3+ 1 ' 1 
z '°-0.7 


3.9 


2.48 


7 


07h 23m 54.14s 


-27°42' 


29.5" 


6.1 


9 9+1.0 


1.2 


1.20 


8 


07h 19m 30.22s 


-27°45' 


42.5" 


4.1 


8-8t 3 2 ; 4 4 


0.5 


0.44 


9 


07h 19m 37.14s 


-27° 33' 


12.0" 


5.2 




3.7 


0.60 


10 


07h 19m 25.72s 


-27°31' 


37.0" 


5.8 


u - M -0.3 


12 


0.36 


11 


07h 19m 18.04s 


-27° 25' 


15.4" 


5.0 


i 7+0.9 


5.9 


0.97 


12 


07h 19m 41.92s 


-27°39' 


58.1" 


5.0 


1 6+ 1 ' 2 
1 '°-0.7 


6.4 


1.02 



TABLE 5 

UVOT photometry for XRT detections in January event 



XRT Source No. 


00031576001 
white 


00031577001 
white 


00031577002 
white 


00031578001 
u 


00031578002 
u 


00031592001 
uvm2 


1 

2 
3 
4 
5 
6 
7 
8 


20.09±0.06 


17.70±0.04 
19.50±0.06 


17.78±0.04 


>21.53 
>21.55 


>21.33 
>21.27 


19.83±0.09 
18.36±0.04 








>21.52 


>21.40 





Note. — Column headings designate UVOT fields and filters used and table entries are observed magnitudes. No XRT sources fell 
within the UVOT field of view for observations 00031575001, 00031575002, 00031576002, 00031579001 or 00031579002. 



TABLE 6 

UVOT photometry for XRT detections in September event 
XRT Source No. 00031825001 00031825003 00031826001 00031826003 



u u white white 

T Z Z ~ Z 

2 >20.73 >21.21 

3 14.11±0.02 14.11±0.02 

4 >20.76 >21.22 

5 19.76±0.17 20.67±0.23 

6 >20.79 21.10±0.33 
7 

8 

9 ... ... 20.75±0.36 20.32±0.12 

10 ... ... 18.76±0.07 18.90±0.04 

11 ... ... >20.79 

12 ... ... >20.84 >21.65 



Note. — Column headings designate UVOT fields and filters used and table entries are observed magnitudes. 



didate event Swift was nominally going to observe five 
0.4° x 0.4° fields. The way the end-to-end search pipeline 
was constructed, only the most probable tiles according 
to the sky- map p m (Q) were observed as Swift fields. This 
resulted in a natural selection bias. As a result, the posi- 
tion correlation term, A cor (f2), had a negligible effect in 
separation of real events from background events. After 
verifying this through simulations we drop it from the 
right hand side of Equation ^ , and the final expression 
for the detection statistic for the joint LIGO- Virgo and 
Swift search becomes 

ftoint(?7, S) = pgw(v) + Pbm(S) (3) 



where pgw{v) an d Pem{S) are logarithms of Aqw(^) and 
Aem('S') respectively. 

We simulated the search by processing a population 
of model GW signals which were paired with plausible 
X-ray fluxes. The set of GW injections was the s ame 
with the one used previously (jAbadie et al.l l2012ah for 
methodological studies of joint GW-EM observations. 
These injections sampled the known galaxies within 50 
Mpc according to the GWGC and were weighted to re- 
flect each galaxy's blue light luminosity. Their intrinsic 
strength (at the source) spanned ad hoc standard candle 
values over 3 orders of magnitude. For the purpose of 
this analysis and to approximate a more realistic distri- 
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TABLE 7 

Potential Variable Stars in UVOT Data 



Variable Source No. 


Right ascension (RA) 


Declination (dec) 


Filter 


Mag 


cm 


Variability 




(J2000) 


(J2000) 










VI 


05:56:22.06 


-40:58:38.3 


white 


19.77 


0.018 


8.03 


V2 


05:56:15.71 


-40:56:06.8 


white 


20.51 


0.026 


9.73 


V3 


05:56:04.51 


-40:51:50.3 


white 


20.50 


0.020 


7.14 


V4 


05:55:54.64 


-40:59:55.3 


white 


19.29 


0.017 


6.35 


Y5 


05:55:43.65 


-40:55:51.9 


white 


21.04 


0.029 


6.88 


V6 


05:55:50.23 


-40:45:40.8 


white 


19.38 


0.063 


12.10 


V7 


05:55:05.25 


-41:04:07.9 


white 


16.53 


0.008 


14.93 


V8 


07:24:20.67 


-27:33:36.5 


u 


17.26 


0.013 


9.67 


V9 


07:23:50.86 


-27:30:24.8 


u 


17.82 


0.024 


8.74 


V10 


07:23:39.28 


-27:25:38.5 


u 


19.42 


0.030 


7.29 


Vll 


07:19:08.99 


-27:37:35.0 


white 


16.27 


0.044 


14.50 



bution of events that are relatively close to our detection 
threshold, we imposed that they fall below detectability, 
which implies r\ < 3.5, at distances outside the 50 Mpc 
range. For the low latency search in 2009-2010, 3.5 was 
the typical threshold value for rj. Translating this value 
into strain at the detector and ultimately at the source 
depends mildly on the waveform morphology and polar- 
ization state of the GW burst and rather stro ngly on the 
frequency content of it (|Abadie et alj l2012bh . In order 
to set the scale, at 50% detection efficiency during the 
2009-2010 run and for GW bursts with energy content 
near 150Hz this value of rj corresponds to an isotropic 
energy at the level of 5.6 x 10 _2 M Q c 2 emitted in GWs 
from a hypothetical s ource at the Virgo cluster (16 Mpc) 
(jAbadie et al.ll2012bO . 

All such simulated events were added to LIGO and 
Virgo data and analyzed as in the actual search. We 
used these simulated GW signals to compute the coher- 
ent network amplitude, "q, and its probability distribu- 
tion, p{rj\signal), which is needed for calculation of the 
joint statistic in Equation ^ and Equation ([3]). Models 
for X-ray counterparts were based on GRB afterglows ob- 
served by Swift. In order to set the scale of possible X-ray 
fluxes for counterparts, we considered sever al short hard 
GRB s and some bright and dim long GRBs (jZhang et al.l 
2009). The host galaxies of all selected GRBs had z < 1. 
We sampled the observed X-ray afterglow light curves for 
the observed GRBs at different time lags (~ 10 4 s, 10 5 
s, 10 6 s) relative to the time of arrival of the burst. For 
a possible X-ray counterpart to GWs from a source at 
50 Mpc away this analysis gave us a wide range of flux 
values, S^oMpc, from 10~ 14 erg s _1 cm -2 to 10 -8 erg s _1 
cm~ 2 . In the absence of any other guidance, we per- 
formed simulations for each order of magnitude in that 
range. For a given S'soMpc every GW signal was paired 
up with a corresponding X-ray flux, which was scaled up 
(as it was positioned anywhere within the 50 Mpc range) 
according to the distance to the source. 

Computation of the joint detection statistic, Equation 
(|3"T). also requires estimates of the background noise in 
GW detectors and flux distribution of serendipitous X- 
ray sources. Background noise in GW detectors is domi- 
nated by high amplitude instrumental artifacts. It is typ- 
ically estimated by time-shifting data from one detector 
with respect to the other. In our simulations we used 
estimates for background noise from the coherent Wave- 
Burst search for GW bursts with LIGO and Virgo in their 
2010 science run (jAbadie et alj2012b| ). In the EM sector, 
on the other hand, serendipitous X-ray sources observed 



in coincidence with a GW signal are the main source 
of false alarms. As in assessing the background in the 
XRT analysis we presented in Se ction 4.1, we used a lso 
here the 2XMMi-DR3 catalogue (| Watson et all 120091) of 
serendipitous sources for the estimation of the flux distri- 
bution for such sources. For both types of backgrounds 
we fitted analytical models to the data and computed 
p(?7|noise) and Pq 1 (S) appearing in the definition of the 
joint likelihood ratio in Equation ((2J. 

For each simulated GW signal complemented by an X- 
ray counterpart of a given flux, the last step of the analy- 
sis involved the calculation of the joint detection statistic 
Pjoint as given by Equation ([3]). Using estimates for the 
GW and X-ray backgrounds, the false alarm probability 
(FAP) of observing a background event with joint statis- 
tic Pj oint > Pjoint in a month long search was computed. 
This defined the significance of the observed event. Fig- 
ure |3] shows the efficiency in detecting these simulated 
GW-EM signals as a function of FAP, Pq, in a joint 
LIGO- Virgo-Swift search using five or ten Swift fields 
and for a wide range of X-ray counterpart fluxes. The 
efficiency is defined as the fraction of simulated signals 
with FAP, Pq < Po- For comparison, we also plot the 
efficiency curve for the GW only search which does not 
use any (Swift) EM follow-up. In the rare-event region 
below FAP of 10 -4 (i.e., corresponding to below ~ 4<r for 
the case of Gaussian statistics) one can see how at fixed 
event significance the efficiency can improve by a signif- 
icant factor depending on the associated EM flux that 
is measured. As expected, the gain increases with the 
brightness of the X-ray counterpart. It reaches a satura- 
tion point at which roughly one-third and one-half, re- 
spectively, of the signals are detected in the searches with 
five and ten fields observed by Swift. This is determined 
by the number of signals for which their true location 
overlapped with the five (or ten) most significant tiles of 
the sky-map. Only these signals were observed by Swift. 
The rest of the signals were "missed" in the follow-up and 
therefore did not benefit from it. Observing more fields 
with Swift increases the chances of locating the X-ray 
counterpart, but at the same time has the negative effect 
of increasing chances of accidental detection of serendip- 
itous X-ray sources (background) . Figure [3] shows that 
the net gain is noteworthy if ten instead of five fields were 
observed by Swift. In general, an X-ray telescope with a 
wide field of view would be optimal for the purpose of the 
joint search. We should note though that for such a tele- 
scope the position correlation term, A cor (r2), in the joint 
likelihood ratio Equation ([2]) may become important and 
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Fig. 3. — Efficiency as a function of false alarm probability for the 
joint LIGO- Virgo and Swift search. The solid (dotted) curves rep- 
resent performance of the joint search with five (ten) fields observed 
by Swift for various values of the flux of an X-ray counterpart (in 
units of erg s _1 cm~ 2 ) at a distance of 50 Mpc, SsoMpc- The 
dashed line is the curve for the GW only search. 

should be included in the analysis. 

6. SUMMARY AND DISCUSSION 

During two periods in late 2009 and 2010 the LIGO- 
Virgo GW interferometric network sent out low latency- 
candidate GW events to partner observatories for rapid 
follow-ups in various EM bands. Two such events were 
followed up as ToOs by Swift. One of the events fol- 
lowed up by Swift was ultimately revealed to be a blind 
injection ar tificially inserted int o the data as a test of 
the system (jAbadie et al.ll2011cfl and the other was a re - 
duced threshold test event (jAbadie et al.lUoilbl I2012bl ). 
Prompt analysis of both the XRT and UVOT data ob- 
tained from the 7 total fields observed showed results 
consistent with expectations for serendipitous sources. 
Given the lack of EM candidates standing out above 
background, these particular observations do not increase 
our confidence in the validity of the GW transients as es- 
tablished by the GW detectors alone. 

Combining GW and EM astronomy will be pivotal 
in maximizing the science in the advanced detector era 
of gravitational interferometers; it may not only in- 
crease our confidence in the detection of GWs but also 
complete our unde rstanding of t he as t rophysics of the 
observed systems (IBloom et al.l 120091: IPhinnev I 120091: 
iStamatikos et al.l 120091: iMetzger fc Berger 1 120121) . Our 
prototype observing program and end-to-end analysis 
has been the first step in joint X-ray and GW obser- 
vations. We demonstrated their feasibility and the con- 
siderable added value joint observations bring. Improve- 
ments to these will continue to be made in the future on 
both the EM and GW side. The relatively narrow FOV 
of instruments such as Swift with respect to the limited 
pointing resolution abilities of GW interferometers makes 
identifying the position of the source on the sky non- 



trivial. A possible fourth detector site in India, Japan 
or elsewhere and continued refinements in source local- 
ization algorithms are likely to reduce the sky-position 
error area. On the EM side, a more highly optimized 
faint source detection scheme for XRT transients might 
yield improvement in EM sensitivity. In late 2011 Swift 
implemented on-board software changes to allow auto- 
matic scheduling sequences of partially overlapping XRT 
FOV exposures in response to ToO observation requests 
for targets with position uncertainties larger than the 
FOV- this will assist the follow-up of GW targets. A 
significant role will also be played by prompt follow-up 
campaigns in the optical band that may provide rapid 
sub-degree source localizations. Such localizations may 
facilitate the subsequent follow-up with narrow FOV in- 
struments including Swift, thus significantly improving 
the chances of capturing the X-ray signatures of GW 
sources. 

In order to carry out multimessenger astrophysics with 
GWs, it will be extremely important to have Swift and/or 
Swift-like satellites capable of rapid pointing, multi- 
wavelength observations and of as wide a field of view as 
possible operating concurrently with the advanced GW 
detector network later in this decade. Maximizing the 
science from GW astronomy will require sensitive partner 
instruments all across the EM spectrum. The successful 
completion of this end-to-end program of EM follow-ups 
by Swift and other observatories during the most recent 
science runs of the LIGO- Virgo network provides confi- 
dence that such joint observations will be both techni- 
cally feasible and scientifically valuable endeavors in the 
future. 
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APPENDIX 

DERIVATION OF THE JOINT LIKELIHOOD RATIO 

In this section we derive from first principles Equation (2) that we used in Section 5 in order to establish the detection 
statistic of the joint GW-EM search. Using Bayes' theorem, the probability skymap for a potential GW source, p m (f2), 
can be written as 

Pm (0 | r,, signal) = p{r] 1 ^^"j Sigiml) (Al) 

p(rj | signal) 

where p(rj | fi, signal) and p{rj | signal) are the conditional probabilities to measure a GW signal with coherent network 
amplitude r\ in the case of a source located at f2 and a source with unknown location, respectively; p(fl | signal) is the 
prior probability distribution for source location, which in this search is determined by the distribution of galaxies in 
the GWGC catalog and distance weighting, see Equation ([1}. 



The joint likelihood ratio is defined as 



c , _ pijh S, Q | signal) _ Jp( V \ Q\ signal)p(S, n | S', SI' , signal)^' | signal)^' | signal) dSl' dS> 

A jointW , - p ^ ^ n ! no . se) - p ^ ! noise)po{s) m 

where Po(S) is the probability density of observing an accidental, serendipitous X-ray source which is not correlated 
with the GW signal, and we introduce p(S, SI | S', SI' , signal), the probability distribution of flux, S, and X-ray source 
location, SI, as measured by Swift for a source whose true flux at the detector and location are S' and SI' , respectively; 
p{S' | signal) is the prior probability distribution for the flux of an X-ray counterpart to the GW signal. Integration 
in the numerator is performed over all possible values of flux and sky locations. Note that, although we assume the 
strength of the GW signal to be uncorrelated with the flux of the X-ray counterpart, we demand that both types of 
signals originate from the same sky position, SI' . This enforces correlation between measured locations of GW and its 
X-ray counterpart. 

Using Equation (|A1|) we can express the joint likelihood ratio in terms of the skymap, p m (Sl), and GW likelihood ratio, 
Agw(^) = p(v I signal) /p(rj | noise): 

J Pm (Sl')p(S, Q | S', SI', signal)?^' | signal) dST dS> 

Aj oint (r),S,n) = AgwW- 77^ • (A3) 

Po(S) 

This expression can be simplified further by postulating p(S, SI 1 S' , SI' , signal) = 5(S—S')8(Sl—Sl') and p(S' | signal) = 1, 

A ioint (r},S,n) = A GW (r?)^r- ( A4 ) 

Defining Aem(<5) = Pq {S) and A cor (f2) = p m (Q), we arrive at the form for the joint likelihood ratio given in Equation 
Ipj). We stress that the particular form of Aem(S) and A cor (J7) is a consequence of the simplifying assumptions 
about Swift's ability to measure the flux and location of an X-ray source that are justified in Section [5] In general, 
these quantities are going to be nontrivial ratios of likelihoods estimating odds of an X-ray candidate source to be a 
counterpart to GW based on its brightness and location. 
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